Low-g dispersion of transverse acoustic phonons in YBa2Cu307_j 
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Abstract 

High-resolution (AE=1.3 meV) inelastic X-ray scattering (IXS) in optimally doped untwinned 
YBa2Cu307_,5 single crystals show that the dispersion of transverse acoustic phonons in the a- 
direction remains linear with no softening due to bending modes in the long wavelength limit, down 
to g=0.02. The sound velocity is found to be 2700 ± 70 m s" 1 , in excellent agreement with ultrasound 
measurements, and unchanged upon cooling from room temperature to 98 K. These results rule out a 
strong coupling between bending dynamics of the Cu02 planes and charge carriers in cuprates. 
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Recent experimental studies renewed the interest in the possible role of lattice dynamics 
in the superconducting state of cuprates. In particular. Inelastic Neutron Q, ^ and X-ray jjj 



Scattering (INS, IXS) studies support early evidence |4| of ubiquitous phonon instabilities in 
the superconducting CuC>2 planes at commensurate in-plane wavevectors q. Further progress in 
the field requires to establish whether these instabilities affect, or merely coexist with, the su- 
perconducting state. On one hand, it has been established that lattice instabilities are common 
in perovskite-like compounds, independent of superconductivity A coexistence between 
these instabilities and superconductivity is supported by a recent study on short-range oxygen 
orderings in YB^CuaOy-^ jjj. The authors found that incipient lattice modulations are in- 
dependent of electronic doping and, hence, of the superconducting critical temperature, T c . In 
La2_ 3; Sr a ;Cu04 , INS data support a picture of incipient stripe order phase with the same peri- 
odicity in superconducting and non-superconducting samples Q. On the other hand, previous 
theoretical studies indicate that finite q anomalies, a signature of incipient lattice modulations, 
may compete fl with superconductivity but also enhance Q the electron-phonon mediated 
pairing, provided the modulations are dynamic and the metallic state preserved. 

In order to clarify this point, in this work we consider a new type of lattice instabilities 
corresponding to bending modes, with uj, q ~ 0. As we shall discuss below, these modes 
are potentially interesting for enhancing T c , whenever a sizeable softening of these modes is 
concomitant to a strong coupling with the charge carriers. In cuprates, as well as in most 
systems, the q ~ region of the phonon dispersions was extensively investigated for the optical 
modes, but remains unexplored for the acoustic ones Q]. To our knowledge, this incompleteness 
of available data is due to the experimental difficulty of combining a high resolution in energy, 
E = hw, with a high resolution in momentum transfer, q. Thus, an adequate experimental set 
up requires simultaneously a highly monochromatic and collimated neutron or X-ray beam, an 
energy analyzer with high energy and momentum resolution and high-quality single crystals 
with a low degree of mosaicity. To explore the above region of phonon dispersion for acoustic 
modes would be important to gain a full insight into the complex dynamics of cuprates, thus 
enabling more stringent tests of competing theories of high-T c superconductivity. 

In this context, the in-plane transverse acoustic (TA) mode is particularly interesting. For 
this mode, a bending dynamics with a softening at sufficiently low q < 0.1 in reduced lattice 
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units (r.l.u.) is expected in layered compounds and would be favored by the covalent in-plane 
Cu-0 bonds of cuprates. According to the standard theory of elasticity in continuous media 
~y\ . assuming isotropic layers, the expected dispersion relation is as follows: 



Utql + Ulql + 1 ql (1) 



where uj is the wave frequency, 7 is the bending stiffness constant of the layer, the suffixes || and 
_L indicate the in- and out of-plane quantities, U = C44/ q would be the square of the velocity 
of sound in the absence of the last term, C44 is the layer shear elastic constant and q is the 
volume density. A bending-induced softening of in-plane TA modes in agreement with the above 



graphite intercalated compounds 



[13], 



expression was experimentally observed in graphite 
and layered transition metal dichalcogenides like 2H-M0S2 |14( . Such peculiar lattice dynamics 
greatly enhances the electron-phonon interaction without invoking any non-linear phenomena. 
Taking into account Eq. (1), in first approximation, the vertex part of the interaction takes 
the following form Q|: 

where w expresses the strength of the electron-phonon coupling, e is an arbitrarily small quan- 
tity and the other symbols are as above. From Eq. (J2J), it is apparent that the vertex tends 
to diverge at small q values for large 7/t/ 2 ratios, i.e. if the mode energy tends to zero, as 
experimentally found in the aforementioned layered compounds where the bending (or shear) 
dynamics dominates. This situation may be relevant to the superconducting state of cuprates, 
since early experimental studies put into evidence a strong coupling between the superconduct- 
ing order parameter and shear lattice distortions in La2_ z Sr x Cu04 and YE^C^C^^ [16]. 
To verify this possibility, in the present work we searched for a possible softening of the in- 
plane TA mode in the hitherto unexplored low-g region. Using IXS with a very high resolution 
in energy, AE—1.3 meV, we measured the dispersion relation from g=0.3 down to g=0.02 in 
YBa2Cu307_5 single crystals. In the q < 0.1 region, the data show no deviation from linearity 
either at room temperature or at 98 K, the latter temperature being well inside the region of su- 
perconducting fluctuations. These results rule out a strong coupling between bending dynamics 
of the CuC>2 layers and high-T c superconductivity in cuprates, and furthermore demonstrate 
the feasibility of accurate measurements of very low-g modes by means of IXS. 



The experiment was carried out on the undulator ID28 beamline at the ESRF. We em- 
ployed beam spot sizes as small as ~ 0.2x0.1 mm 2 and two different Si (11 11 11) and (13 13 
13) monochromator Bragg reflections Using the ID28 6m-long arm spectrometer, these 

reflections provided an overall energy resolution of 1.7 and 1.3 meV, respectively 19]. We in- 
vestigated the in-plane TA dispersion along A (a or [£,0,0] direction) with out-of-plane ([0,0,£]) 
polarization. One can verify that the q resolution for this phonon, AQ a , is limited by two 
independent factors: 1) the mosaic spread of the crystal, A$ c , convoluted with the divergence, 
AKi, of the incident beam, the convolution being measured experimentally by the rocking 
curve width, Au; 2) the transverse component, A2# t , of the width, A28, of the scattering solid 
angle. In a typical INS experiment, AQ is usually dominated by the former factor. In the 
present IXS experiment, one can significantly enhance AQ a thanks to the very small beam 
divergence (120 /irad horizontal x 60 /irad vertical) characteristic of the synchrotron undulator 
sources and optimized beamline optics, provided a high quality crystal with a sufficiently low 
degree of mosaicity is used. This is the case of the present experiment. We measured two 
YBa2Cu307_,5 single crystals grown by the flux method, as described elsewhere j^. The 
crystals are small platelets of about ~ 1 x 1 x 0.1 mm 3 size selected from two batches pre- 
pared under nominally the same conditions. Oxygen post-annealing at 500 °C for two weeks 
enabled to achieve a full-oxidation (5 f=s 0.1). The degree of mosaicity was very low in both 
crystals, as indicated by the small full widths at half maximum of the (006) rocking curves, 
Au;=0. 009-0. 018° for both samples, that are close to the horizontal beam divergence, 0.007°. 
This gives a mosaic spread of 0.005-0.016°. Using the Si (13 13 13) reflection, with A.;=0.4824 
A, it turns out that the contribution of Auj to AQ is in the 0.0003-0.0006 r.l.u. range in the 
(0,0,6) Brillouin zone. We used two configurations of horizontal slits with horizontal openings of 
0.190 and 0.095°, the latter for the measurements at q points near the zone center. To minimize 
the spurious contribution from the (0,1,0) direction near the zone center, the opening of the 
vertical slits of the analyzer was reduced to decrease the sagittal scattering angle from 0.57 to 
0.19°. A straightforward calculation shows that, under the above experimental conditions, AQ 
along the longitudinal, [£,0,0], and transverse, [0,0, Q, directions in the (0,0,6) Brillouin zone 
is in the 0.001-0.002 and 0.04-0.08 r.l.u. range respectively. Hence, in the present experiment, 
AQ is dominated by the divergence of the scattered beam, i.e. by the analyzer slit opening. 
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FIG. 1: Energy scans near (0,0,6) for various values of the in-plane momentum transfer £ along the 
a-axis (indicated in the legend) for one of the two YE^CusO^-^ single crystals. These data were 
obtained using the Si (13 13 13) reflection, which provides an overall energy resolution of 1.3 meV. 



In the case of the Si (11 11 11) reflection, with Aj = 0.5701 A, the estimated AQ values are 
increased by an amount comparable with the 16% increase of Aj of this reflection with respect 
to the (13 13 13) one. 

The phonon dispersion was measured in a fixed-Q configuration at various (£,0,6) points 
of the Brillouin zone in the vicinity of the (006) Bragg reflection. The latter reflection is 
the most intense one of the (00/) type in fully oxygenated YBa2Cu307_,5 , thus providing the 
largest photon-phonon scattering cross-section for low-g acoustic phonons. The same runs of 
measurements were carried out on the two crystals. In Fig. ^ a series of room temperature 
energy scans at various values of phonon momentum, £, is reported. For all £ values, one clearly 
distinguishes the triple peak structure consisting of the central peak, that corresponds to diffuse 
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FIG. 2: Comparison between the TA phonon peaks at 300 and 98 K measured using the Si (11 11 
11) and (13 13 13) reflections. Solid lines are a fit of the data as described in the text. 

elastic scattering, and of the two lateral phonon, Stokes and anti-Stokes, peaks. The Si (11 11 
11) data enabled to resolve the phonon peaks down to £ = 0.04. At lower £ values, these peaks 
appeared as too small shoulders of the central peak and could not be resolved. To resolve them, 
it was necessary to use the (13 13 13) reflection and closer analyzer slits, as described above. 
The relatively low intensity of the central peak as compared to the phonon peaks indicates 
that the amount of disorder, responsible for the diffuse scattering, is limited in the crystal. For 
every energy scan, the above triple peak structure was fitted using the standard constraints 
regarding the relative position and intensity of the two lateral Stokes and anti-Stokes peaks 
set by the Bose-Einstein statistics. The phonons peaks and the diffuse elastic scattering peak 
were modelled by harmonic oscillators and a delta function respectively. The sum of these 
contributions was convoluted with the instrumental resolution function as described in Refs. 
1^, Four examples of fitted energy scans are shown in Fig. |21 The agreement within the 
error bars between the experimental data and the model is evident. The phonon energy values 
extracted from all scans yield in a straightforward manner the dispersion reported in Fig. EJ 
where the results obtained at room temperature and 98 K are summarized. One notes that the 
data of the two samples are reproducible within the experimental uncertainty The accurate 
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IXS measurements enable us to firmly conclude that the dispersion is linear down to g=0.02. A 
linear regression in the q < 0.1 range of the data of sample 810R, which gave the most accurate 
results thanks to the lower mosaicity, yields a transverse sound velocity along the a direction of 
2700 ± 70 m s _1 . This value is in excellent agreement with the value of 2710 m s _1 derived from 
a previous ultrasound measurement of C55 |22J, the relevant elastic constant for orthorhombic 
YBa2Cu307_5 . Good agreement is also found with the value of 2620 m s _1 obtained from a 
shell model calculation of C55 j^ . The discrepancy in the latter case is compatible with the 
approximations of the model. We finally fitted the same data set using a sine function to take 
into account the levelling off of the mode frequency at the zone boundary. This fit yields a 
slightly higher velocity, 2810 m s -1 , which suggests that the linear fit may underestimate the 
slope of the dispersion curve even at q < 0.1. In any case, it is apparent from Fig. |3]that the 
present measurements agree with the ultrasound ones within the error bars. 

Neither a hardening nor a softening of the mode is found upon cooling from room temper- 
ature down to 98 K, as it is noted from Fig. Eland, in more detail, Fig. |2j It is recalled that 
the latter temperature is well inside the region of thermodynamic fluctuations of the supercon- 
ducting order parameter in YBa2Cu30"7_5 . It follows that no change of the lattice dynamics 
regarding the in-plane TA mode occurs in presence of the superconducting instability. By as- 
suming the relevance of the electron-phonon interaction to the pairing mechanism in cuprates, 
one concludes that this mode would not be strongly involved in the pairing process. At best, 



according to a well established result of the BCS theory |24(, a renormalization of the phonon 
frequency and lifetime is expected at lower temperatures as a general consequence - rather than 
a cause - of the superconducting transition that affects all phonon modes. 

In conclusion, using IXS with very high resolution in energy and momentum, we performed 
the first measurement of in-plane TA phonon dispersion in the very low g-limit on cuprates. 
Stringent conditions regarding the monochromaticity and the collimation of the beam combined 
with the low degree of mosaicity of the samples enabled us to measure reproducibly the dis- 
persion along the a*-axis in two untwinned YBa2Cu30"7_,5 single crystals. Upon cooling from 
room temperature to 98 K, well inside the region of superconducting fluctuations, no anomaly 
was found down to q = 0.02, the dispersion remains linear within the experimental uncertainty 
and the corresponding sound velocity is unchanged. Its value is in agreement with previously 
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FIG. 3: Dispersion relation of the in-plane TA phonon at room temperature and 98 K for the two 
YBa2Cu30"7_,5 crystals labelled as 810R and 811R. The right panel shows in detail the very low q 
data for the former. Vertical lines indicate the energy and momentum resolutions that depend upon 
the Si (11 11 11) or (13 13 13) reflection used, slit openings and crystal mosaicity (see text). 

reported ultrasound measurements and shell model calculations. The absence of any anomalies 
or temperature-induced mode hardening or softening in the region of superconducting fluctu- 
ations indicates that, in cuprates, the bending dynamics and the coupling between in-plane 
TA acoustic modes and the electronic degrees of freedom are negligible. This result should set 
stringent conditions for any phonon-based model of high-T c superconductivity in these materi- 
als. In addition, our experiment shows the suitability of IXS to accurately measure low-energy 
excitations in the q=0 limit, provided single crystal with a low degree of mosaicity are used. 

This work has been supported by the European Synchrotron Radiation Facility (Experiment 
HS-2445). The authors acknowledge the technical assistance of T. Besagni and P. Ferro. 
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